Hadron-nucleus scattering in the local reggeon model 
with pomeron loops for realistic nuclei 



M.A.Braun, A.Tarasov 
S.Peterburg State University, Russia 

March 26, 2010 



Abstract 

Contribution of simplest loops for hadron-nucleus scattering cross-sections is studied in the 
Local Reggeon Field Theory with a supercritical pomeron. It is shown that inside the nucleus the 
supercritical pomeron transforms into a subcritical one, so that perturbative treatment becomes 
possible. The pomeron intercept becomes complex, which leads to oscillations in the cross-sections. 

1 Introduction 

High-energy hadron-nucleus scattering has long been studied in the framework of the local Regge- 
Gribov theory with a self-interacting supercritical pomeron. The first successful description of the 
total cross-sectuion in this approach was presented by A.Schwimmer in 1975 [1] who performed 
summation of fan diagrams in the limit when the pomeron slope a' was zero. With the advent of 
the QCD attention has been shifted to self-interacting non-local BFKL pomerons. In this approach 
an analog of the Schwimmer model has been constructed in the form of the Balitsky-Kovchegov 
evolution equation [2, 3, 4]. However QCD can reliably describe only the hard region of the 
dynamics. Soft processes, contributing to the bulk of the total cross-section, are better described 
by the old-fashioned local reggeon theory, so that the latter remains a useful tool in spite of its 
venerable age. 

Both the old local and new non-local pomeron models have been fully considered and applied 
only in the quasi-classical tree approximation. Loops have been neglected in both formalisms. This 
approximation can be justified if the parameter 7 = Aexp Ay is small, with y the rapidity and A 
and A the pomeron intercept and triple pomeron coupling. Then for a large nuclear target, such that 
^1/3^ the tree diagrams indeed give the dominant contribution and loops may be dropped. 
However with the growth of y the loop contribution becomes not small and this approximation 
breaks down. 

Full calculation of the loop contribution seems to be a formidable task for the non-local QCD 
pomeron. So it seems to be worthwhile to start with the Local Reggeon Field Theory (LRFT) 
with a supercritical pomeron. Such a study, apart from its possible lessons for the modern QCD 
approach, also has an independent value, since, as mentioned, the old LRFT with phenomenological 
parameters describes the soft dynamics of high-energy strong interactions not so badly. Much work 
has been done on the influence of loop diagrams in the zero-dimensional LRFT, both theoretically 
old ago in [5, 6, 7, 8] and numerically recently in [9]. This influence turns out to be decisive for 
the asymptotic behaviour at large energies and transforms the supercritical pomeron into a weakly 
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subcritical one with the effective intercept oc — exp(l/A^) where A is the triple-pomeron coupling 

assumed to be small. 

Unfortunately generalization of these beautiful results to the realistic case of two transverse 
dimensions is prohibitively difficult. To start with, one is forced to introduce a non-zero value of 
the slope: otherwise the loop contribution is divergent in the impact parameter. However even in 
the tree approximation solution of the model with a' 7^ is only possible numerically. Second, the 
model in = 2 needs rcnormalization in the ultraviolet. And, most important, the method used 
to solve the model in dr = 0, which is to study the corresponding quantum-mechanical system and 
the relevant Schroedinger equation, is inapplicable in the realistic case, since instead of ordinary 
differential equations one arrives at equations with functional derivatives. In fact summation of all 
loop contributions is equivalent to a complete solution of the corresponding quantum field theory, 
the task which seems to be beyond our present possibilities. So at most one can hope to obtain 
some partial results which might shed light onto the properties of the model with loops. There were 
several attempts to study the high-energy behaviour of the LRFT with a supercritical pomeron 
using different approximate techniques and giving contradicting results [10], [11]. 

In our previous study, instead of trying to solve the model for the purely hadronic scattering we 
considered the hadron-nucleus scattering and propagation of the pomeron inside the heavy nucleus 
target [12]. Moreover to avoid using numerical solution of the tree diagrams contribution with 
diffusion in the impact parameter, we concentrated on the case of a constant nuclear density which 
allowed to start with the known analytical solutions. We have found that the nuclear surrounding 
transforms the pomeron from the supercritical one with intercept e > to a subcritical one with 
the intercept — e. Then Regge cuts, corresponding to loop diagrams, start at branch points located 
to the left of the pomeron pole and their contribution is subdominant at high energies. As a result 
the theory aquired the properties similar to the standard LRFT with a subcritical pomeron and 
allows for application of the perturbation theory. 

However the adopted approximation of a constant nuclear density in the whole space was ob- 
viously too crude. In fact it did not allow to make convincing conclusions about the physical 
cross-sections, so that we had to recur to some weakly supported guesses on this point. In this 
paper we try to improve our treatment and consider realistic nuclear targets which occupy finite 
space volume. This enables us to predict the influence of loops on the total cross-section and study 
applicability of the perturbative approach to hadron-nucleus scattering with pomeron loops. 

Our results depend on the behaviour of the nuclear profile function T{b) at large values of the 
impact parameter b. For a finite nucleus with T{b) = at 6 > Ra all our conclusions found in 
[12] for constant T{b) remain valid and one can apply perturbative methods to study the loop 
contribution. However for an infinite nucleus, for which T{b) exponentially falls with b but never 
vanishes at finite b, the loop influence is perturbative only inside the nucleus, at not too large b. 
At very large b, well outside the bulk of the nuclear matter, the loop contribution blows up and 
the cross-section becomes unperturbative. 

The paper is organized as follows. In the next section we introduce the model and reformu- 
late it to describe the loops in the nuclear surrounding. In Section 3 we introduce a quasi-local 
approximation to constructively treat the case of realistic nuclei. In Section 4 we formulate the 
appropriate Dyson equation for the amplitude with loops. Our numerical results for total pPb and 
pCu cross-sections are presented in Section 5. Finally Section 6 draws some conclusions. 
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2 Local pomeron in the nuclear field 

The LRFT model is based on two pomeron fields (j){y, b) and (/>^(y, h) depending on the rapidity y 
and impact parameter 6, with a Lagrangian density 

L = Lo + A</.t</.(0 + 0t)+5P</,. (1) 
Here the free Lagrangian density is 

^0 = 0^ 5 -a'Vl + e) = </.t5</. (2) 

where e is the intercept minus unity and a' is the slope. The source term describing interaction 
with the nuclear target at low energies is 

gp{y,h)=gAT{h)5{y). (3) 

where g is the pomeron-nucleon coupling constant and T{h) the profile function of the nucleus. For 
a supercritical pomeron e > and A < 0. 
The classical equation of motion are 

— = -dy(l)^ - a'Vlcj)^ + e<^t + A<^t^ + 2A<^<^t + gp = {) (4) 
o(p 

and 

A T 

= dyc/)^ - a'Vl<p + £(/) + A(/)2 + 2A(/>t(/> = 0. (5) 



Prom the latter equation we find = and the equation for (/>^ takes the form 

dy(f)^ = -a'Vl(f>^ + e<^t + X(f>^\ (6) 

with an initial condition 

4>\y = 0)=gp (7) 

Equation (6) describes evolution of the pomeron field in rapidity and its diffusion in the impact 
parameter inside the nucleus. We denote the solution of the classical equation of motion (6) with 
the initial condition (7) as ^{y, b) 

To go beyond the classical approximation and thus study loops we make a shift in the quantum 
field <^t: 

<^t(y,6) = 0t(y,6)+e(y,6) (8) 

and reinterprete our theory in terms of fields (f) and ^J. In the Lagrangian terms linear in (f) vanish 
due to the equation of motion for ^ and we obtain 

L = 4iS + 2X0cf> + A^.^2 + ^^t^^^t + (9) 

This Lagrangian corresponds to a theory in the vacuum with the pomeron propagator in the external 
field f{b,y) = 2\ay,b) 

P = -{S + 2X0-\ (10) 

the standard triple interaction and extra interaction described by the term A^(/>^. This new inter- 
action corresponds to transition of a pair of pomerons into the vacuum at point {y, b) with a vertex 
A^(y,6), see Fig. 1. 



3 



Figure 1: The new vertex for twopomeron annihilation, which appears after the shift in field (f)^ 




Figure 2: Some simple loop diagrams for the pomeron Green function 



Loops can be formed both by the standard interaction and the new one. In the latter case they 
are to be accompanied by at least a pair of of standard interactions. Diagrams with a few simple 
loops in the Green function are illustrated in Fig. 2. One immediately observes that a loop formed 
by the standard interaction has the order }? /a' and requires renormalization. A loop formed by 
the new interaction has the order A^/a' and is finite. 

The amplitude is obtained as a tadpole g < (j)i{y,b) >. The simplest diagrams for it contain 
one loop and are shown in Fig. 3 a, b. Diagrams with more loops are shown in Figs. 3 c, d. 

Propagator (10) in the external field f{y,b) = 2X^{y,b) where ^{y,b) is the classical field 0^ 
satisfies the equation 



dPiy,b\y'b') 
dy 

with the boundary conditions 



(e - a'Vl)P{y, b\y' , b') + /(y, b)P{y, b\y' , b') 



P{y,b\y',b') = 0, y-y'<0, P{y' ,b\y' ,b') = d\b - b') 



(11) 



(12) 



In the general case propagator P can only be calculated numerically, just as the external field 
f{y, b). Its analytic form can be found in two cases [12]. 
If the slope a' = then one finds 



f{y,b) 



ea{b)e' 
p{y, b) 



(13) 
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Figure 3: Diagrams with one loop (a, b) and two loops (c and d) for the scattering amplitude 



and 



where 



and 



P(y, b\y', b') = 5\b- b')Poiy, y', b), Po{y, y', b) = e^(f-f')^|^, 



a{b) 



\gp{b) 



> 



p{y,h) = l + a{b)(e'y -I). 



(14) 



(15) 



(16) 



Propagator PQ{y,y',b) corresponds to the zero-dimensional theory with a{b) as an external param- 
eter. It is remarkable that at large y propagator Po{y, y', b) ~ exp(— ej/), that is behaves as the free 
propagator with the opposite sign of e and so vanishes at y ^ oo. In contrast to the free propagator 
it corresponds to a subcritical pomeron. 

The second case which admits an analytic solution for the propagator is that of the nuclear 
matter, that is the case when the profile function is constant in all the transverse space 



Tib) = To. 



(17) 



Then the classical field ^ and the external field / = 2A^ become 6-independent and the propagator 
in this field takes the form 



P{y,b\y',b') 



1 



A-Ka'{y — y') 



y') 4oc'(y-y') 



(18) 



where p{y) is given by (16) with a constant a. Note an especially simple case when a = 1 and 
p{y) = expey. Then the propagator in the external field coincides with the free propagator with 
the opposite sign of e, so that the theory formally corresponds to a subcritical pomeron model with 
an additional interaction shown in Fig. 1 



3 Realistic nucleus 

The full treatment of the scattering on a realistic nucleus with the profile function T{b) which varies 
with b and vanishes as & — > oo requires solution of the equations (6) for the pomeron field and (11) 
for the propagator in the found external field f{y, b). Both tasks can be done only numerically and 
look hardly feasible. 
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To facilitate the problem we shall use the approximation of a weakly varying T{b) . Equivalently 
we assume that the slope a' , which governs diffusion in the impact parameter, is small. Accordingly 
we shall put a' to zero wherever possible. In particular we shall use expression (13) for the field 
at fixed value b of the impact parameter. For the propagator we shall use expression (14) so long 
as it stands alone. However we cannot do this in pomeron loops, where it is squared leading to a 
divergence at a' — t- 0. 

For loops we shall restrict ourselves with the leading term as a' and so use the leading 
approximation for the propagator inside the loop in this limit but not the limit itself. This leading 
approximation obviously coincides with the approximation of a slowly varying nuclear field and is 

given by Eq. (18) with p depending on b in accordance with definition (16). One can see it from 
the following reasoning. Let us seek the propagator satisfying Eq. (11) as a product 

Piy,b\y',b') = Po{y,y',b)P,{y,b\y'b') (19) 

with the known Po{y, y', b) defined in Eq. (14). Then we find an equation for Pi 

Poiy, y', 6)^^i^^^^ = a'VlPoiv, y', b)P,{y, b\y', b). 

Assuming the smallness of a' and non-singularity of Pq as a function of b or equivalently a weak 
dependence of Pq on b we can neglect action of the derivatives in b on Po{y,y',b) to obtain an 
equation for Pi 

dPi{y,b\y',b') . , 
= a V6Pi(y,6|y ,6). 

Solution of this equation brings us to propagator (18). 

With this approximation we now calculate the simplest loops which appear in the theory paying 

special attention to their 6-dependence. 

The simple loop Fig. 2a is given by the expression 

Si(y, b\y', b') = -2X^P\y, b\y', b') + AeS\b - b')5{y - y'). (20) 

The second term comes from the intercept renormalization. Using (18) we find 

Si(,, 61,', 60 = -2A^^ ( + A.^^(6 - 6')<5(y - ,')• (21) 

The exponential factor depending on (6 — 6')^ obviously becomes proportional to 5^{b — b') in the 
Umit a' ->■ 

_J_,-5S.„,^„^J^(6). (22) 

Thus we obtain at small a' 

Si(y, b\y\ b') = -5\b - b') f e'<y-y')P^ + Aed'{b - b')5{y - y'). (23) 

4Tra'{y-y') P^[y,b) 

Note that at fixed b and rapidities Si is finite. Ultraviolet divergence appears in the process of 
integration over rapidity at y — y' = 0. Regularizing this integration by cutting it at y — y' = ymin 
one can express 

Ae = - - — - \n{cRymin) (24) 
47rQ:' 
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with cr a renormalization constant taking place of Ae. 

Now we pass to the simplest new self-mass of Fig. 2b. It does not need renormalization and is 
given by 

S2(y, b\y', b') = 4A3p(y, h\y'b') f dzd^c^z, c)P{y, b\z, c)P{y', b'\z, c). (25) 



With the approximate form (18) for the propagator we find a product of exponential factors de- 
pending on impact parameters 



1 (b b')2 (j, )2 (^/ )2 

e 4a'(y-y') 4a'(j/-2:) Aa'(y'-z) 



{A-Ka'Yiy - y'){y - z){y' - z) 
to be integrated over z and c. Integration over c converts this into 

1 (i,-6')2 (b-b'f I 

{A'Ka'f{y-y'){y + y'-2z) Sttq' ^ ' 

So in our approximation we find the new self-mass 

S,fe,6|,',6') = -S^(6-6')«fe-.'):^^ r^e<->-'>'^ (26) 

Sttq:' Jo y-z p {y,b) 

4 Amplitude at fixed impact parameter 

In the adopted approximation all integrations over intermediate impact parameters are performed 
with the help of 5 functions and at fixed impact parameter b the final amplitude r(y, 6) is a function 
of b presented as a mutiple integral over intermediate rapidities. So the whole picture is local in b. 

The propagator (14) in the external field has the crucial property that it decreases with y as 
y oo: 

P{y, y', fe)(,-,')^oo ^ e-<y-y'^f{y', b). (27) 

Because of this property internal integrations over rapidities do not lead to contributions which 
grow faster than the lowest order term. This ensures convergence of the perturbative expansion at 
small enough values of A 

Introduction of loops into the amplitude A is achieved by the Dyson equation 

A{y,b)=A^''\y,b)+ dyi dy2P{y,yi,b)J:{yi,y2,b)Aiy2,b), (28) 
Jo Jo 

where the lowest order term is 

with p{y, b) given by (16) and the pomeron self mass S(y, b\y'b') presented in our approximation as 

Eiy,b\y',b') = S\b-b')Eiy,y',b) (30) 

In accordance with the convergence of the perturbative expansion, in the first approximation 
we take 

S(2/, b\y', b') = Si(y, b\y' , b') + ^2{y, b\y', b') (31) 
where Si and 112 are the 2nd and 3d order self-mass contributions studied above. We present 

^(y,5)=^(0)(y,5)r(y,5). (32) 
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The equation for r{y) takes the form 

riy,b) = l + Xiiy,b)+X2{y,b), (33) 
where Xi and X2 are parts coming from Ei and S2 respectively. ExpHcitly 

where 

Ana' 

and 



.3 

The part ^2(2/) is 



^, b) = T dyi^i|^^r(yi - z,b). (36) 



X2{y) = --^P^ dyi^^J -^^pHy3,b)e-^'' dy2-^My2,b). (37) 
p{y,b) Jo P^{yi,b) Jo yi-y3 Jyz p{y2,b) 

Eq. (33) contains the renormahzation constant cr as a parameter. To fix it we use the pre- 
scription proposed in [12]. The contribution to A{y,b) of a single loop is 

^«(y,6)=^W(y,6)r«(y,6) (38) 

where 

r«(y,6) = [X,{y,b)+X2iy,b)) (39) 



The right-hand side of Eq. (39) can be calculated and shown to have a finite limit at y — >^ 00 

^«fe.*W = -i|7;(l"f + (40) 



where C^; is the Eiler constant. Remarkably this limit does not depend on b and is universal. So 
one can use it to fix the renormahzation constant by requirement that this limit is zero: 

CR = 2ee^^-^ (41) 

This means that at large rapidities the loop contribution vanishes, so that one defines the pomeron 
intercept from the asymptotic of the amplitude requiring it to coincide with the lowest order term. 

With all this, one should be careful in applying these results for different values of impact 
parameter b. The decrease of the propagator P{y, y', b) obviously starts when p{y, b) begins to rise, 
which requires a{b)e'^y >> 1 or 

a{b) 

It follows that for very small a{b) the actual decrease will start at very large rapidities. At 
a{b)e''y ^ 1 the propagator will have order l/a(6). In other words, at small a{b) we expect the ratio 
r{y,b) to rise with rapidity from unity to a large value l/a(6) at ey = yo ^ — ln(a(6)) and then 
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drop to zero at y >> yo- At its maximum the loop contribution will be greater than the lowest 
order by factor 

^ " 47ra'ea(6) ' ^^^^ 

where e in the denominator comes from the essential interval of integration over intermediate 
rapidities. If the triple pomeron coupling constnt A is small enough this does not prevent using 
perturbation theory since A^/a(6) ~ A. However it docs bring difficulties if one fixes A and then 
varies a{b) to very small values at large b. This point is essential in the application of the model to 
realistic nuclei studied in the next sections. 



5 Numerical results 

The total hadron-nucleus cross-section is given by twice the imaginary part of the hA amplitude, 
integrated over all values of the impact parameter. To take into account unitarity at fixed b we 
glauberize the contribition from a single fan with loops, which implies that the projectile hadron 
may emit any number of such structures. 

a{y) = 2 J d^(l - e3'^(2''^)) . (43) 

To calculate A{y, b) we have solved the Dyson equation (28) for proton-lead and proton-copper 
collisions {A = 207 and 64) at each value of b with the pomeron self-mass given by the sum of 
lowest order loops Si 4- S2. 

We used more or less standard values for the pomeron parameters extracted from the experi- 
mental data on proton-proton scattering 

e = 0.08, a' = 0.2 GeV'^, X = -0A8GeV-^, 5 = 5.94 GeV'^. (44) 

Note that e and g can be rather reliably taken directly from the proton-proton inelastic cross 
seections. The rest of the parameters is not so well established and model dependent. Our set 
agrees with the parametrization used in [13]. With the set (44) the value of parameter ^ which 
controls the relative order of loop contribution and so applicability of perturbative treatment is of 
the order l/a{b). For heavy nuclei the maximal value of a(b) is a(0) ~ 1. So strictly speaking for 
realistic nuclear scattering the used values of A do not lie within the applicability of perturbative 
treatment. This circumstance has to be taken into account in analyzing our results with A given 
in (44), which should be taken as qualitative. 

We considered two models for the spatial structure of the nuclear target. The simplest one 
assumes the nucleus to be a sphere with a constant nuclear density of radius 

Ra = A^/^Ro, ro = 1.15 fm. 

This is the nuclear model most suitable for our approach, since it excludes the region of very small 
nuclear densities ouside the nucleus in which the perturbative treatment is impossible due to large 
values of (. In this case the found hA cross-sections are shown in Fig. 4 as a function of rapidity 
y. To see the influence of loops we also presented the cross-sections without loop contribution, 
which correspond to the monotonous curves. To clearly see the ^-dependence we rescaled the pCu 
cross-section to be equal to the pPb one at small y. As one observes, the loop contribution is quite 
noticeable and leads to oscillations in the cross-sections. The origin of oscillations is understandable 



9 




Figure 4: pPb and pCu cross-sections on the target nucleus as a sphere with a constant density. 
The oscillating (monotonous) curves corresponds to presence (absence) of pomeron loops. The pCu 
cross-sections are multiplied by 2.55 to coincide with the pPb one at small y 



from the fact that with loops the original pomeron pole splits into two complex conjugated poles, 
which remain on the physical sheet due to the absorptive nature of the loop (negative signs in 
(23) and (26)). The A dependence is quite strong and also understandable. With the growth of A 
the damping effect of the nuclear surrounding grows, so that the cross-sections generally become 
relatively smaller, both with and without loops. 

As a second model for the nucleus density we considered the standard Woods-Saxon one. Such 
a nucleus has no finite dimension and we have to vary h from zero to infinity. With the growth 
of h outside the nucleus parameter a(6) becomes exponentially small and formal application of our 
formulas becomes impossible because of the growth of the perturbative parameter C,. To overcome 
this difficulty we take into account that at b outside the nucleus the projectile practically interacts 
with a free nucleon and so its aplitude is wholly determined by the single pomeron exchange without 
any loop insertions, which are assumed to be incorporated into the effective pomeron parameters. 
So we have introduced a cutoff into the loop contributions multiplying them by 

ao + a{h) 

and choosing ao = a(i?^ + A6) with h = 1.5 fm. The resulting pPb cross-sections are presented in 
Fig. 5 again together with the cross-sections without loops, which correspond to the monotonous 
curve. The found cross-sections show the same qualitative behaviour as with the finite nucleus, 
although the absolute values and the amplitude of oscillations are somewhat greater. The cross- 
sections without loops naturally grow faster that for a finite nucleus because of the contribution of 
the nuclear halo. 

The amplitude of oscillations is determined by the value of A. With smaller A's the oscillations 
die out. To see that we repeated our calculations with a smaller value of the triple pomeron coupling 
constant found in [14] A = —0.069, that is nearly ten times smaller than in (44). With such A 
the perturbation parameter drops to C, ~ 1/5 inside the nucleus, so that perturbative expansion 
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Figure 5: pPb cross-sections on the target nucleus with the Woods-Saxon density. The oscihating 
(monotonous) curve corresponds to presence (absence) of pomeron loops 

becomes reasonable. Our results for pPb cross-sections with A = —0.069 are shown in Fig. 6. We 
observe that the influence of loops is quite small. It grows with rapidity and becomes noticeable at 

y ~ 30 -7- 40 only to vanish in the end at still higher rapidities. Comparing with Figs. 4 and 5 one 
also observes that with smaller A the cross-sections rise faster with y due to the behaviour of the 
zero-order term (29). At very small A the damping exersized by fanning is small and A^^^ grows 
like the pomeron itself. 

6 Conclusions 

We studied the contribution of the two simplest loops in the local reggeon field theory with a 
supercritical pomeron in the nucleus. We applied a quisi-local approximation which assumes that 
the slope a' is small and can be put to zero in all places except in the loops where the leading 
term at a' — > has been retained. Our results confirm the conclusion reached in [12] that the 
nuclear surrounding transforms the supercritical pomeron with the intercept q;(0) — 1 = e > 
into the subcritical one with the intercept — e. As a result, at high energies the pomeron Green 
function vanishes and contributions from multipomeron exchanges vanish still faster according to 
the standard predictions for the subcritical pomeron. However this phenomenon becomes effective 
at different energies depending on the nuclear density. At small densities the Green function begins 
to decrease at very high rapidities far beyond our present experimental possibilities. And in the 
intermediate region loop contributions may become large if the triple pomeron coupling constant 
A is not small enough. 

So while the general conclusion about the validity of the perturbative treatment is true, its 
practical realization depends on the value of A and is determined by the perturbative parameter C, 
defined in (42). Unfortunately with the standard values (44) for A parameter C, ^ 1 and perturbative 
treatment is dubious. Higher order loops have to be included into the picture to have quantitative 
results. With smaller values of A the loop contribution becomes small and perturbation approach 
is justified. 
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Figure 6: pPb cross-sections on the target nucleus with the Woods-Saxon density and A = —0.069 
The lower (upper) curve on the right side corresponds to presence (absence) of pomeron loops 

The influence of loops shows itself in oscillations in rapidity of the cross-sections. These oscil- 
lations are due to the complex character of the pomeron pole once loops are included. Due to the 

wrong sign of the pomeron self-mass the pomeron poles does not move to the second sheet of the 
energy plane as with normal particles but stays on the physical sheet acquiring a nonzero imaginary 
part. Of course with very small A oscillations are small and have a very long pariod in y. However 
with greater A they are clearly visible, as illustrated in Figs. 4 and 5 

On the practical side our results point to yet another possibility to experimentally measure the 
value of the triplc-pomeron coupling constant A on which the behaviour of the cross-sections with 
energy critically depends. To exclude the region of impact parameters b well outside the nucleus, 
where our treatment becomes invalid, one may select events in which the collisions occur inside the 
nucleus by requiring the multiplicity of produced particles to be large enough. 

From the theoretical point of view we consider our results to be promising as a basis for treating 
loop contributions in the perturbative QCD. It seems to be advantageous to study them in the 
nuclear surrounding, which makes the perturbative approach at high energies much more tractable. 
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